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Abstract 


A scattering formalism is developed in a multiple scattering model to describe inclusive 
momentum distributions for high-energy projectiles. The effects of final state interactions on 
response functions and momentum distributions are investigated. Calculations for high-energy 
protons that include shell model response functions are compared with experiments. 

Introduction 

A realistic description of galactic cosmic ray transport through bulk shielding requires an 
extensive data base of nuclear interaction cross sections. (See ref. 1.) Traditionally, velocity- 
conserving interactions have been assumed for all ions with a mass number greater than one for 
easier numerical computations. (See ref. 1.) As transport algorithms become more accurate, this 
assumption will likely be removed; the result will be an increased need for nuclear data bases 
that include secondary particle spectra. 

Previously (refs. 2 and 3), we developed a multiple scattering series to describe the energy 
loss spectra of fast ions in nuclear collisions. This multiple scattering series, which describes the 
quasi-elastic peak, was effectively summed by expressing the many-body response functions of the 
target as a convolution of the one-body response function through an energy shift approximation. 
For composite projectiles, incoherent corrections were shown to slightly reduce the cross section 
(refs. 2 and 3), when uncorrelated form factors were assumed. Multiple scattering effects 
were also shown to shift the position of the quasi-elastic peak considerably, which agrees with 
experiments. In references 2 and 3, a two-dimensional representation of the target response 
function was used in the eikonal formalism. In this paper we treat the longitudinal aspects of 
the response; however, we used the eikonal approximation to evaluate the cross section. In this 
paper we also focus on the final state interaction (FSI) between knocked-out target nucleons and 
the target recoil. A treatment of the FSI in inclusive scattering is important for understanding 
scaling phenomena (refs. 4 -6) and transparency (refs. 7-9) in high-energy collisions as well as in 
the intranuclear cascade. 

In inclusive scattering with fast ions, the projectile has left the scattering region before 
the FSI; thus, a simplified treatment of ion effects on the projectile wave function may be 
valid. Focusing on energy losses above low-lying collective states and below pion production 
thresholds, we consider ejectiles of the target using the approach of Horikawa et al. (ref. 10) 
for decomposition of the response function into elastic and inelastic FSFs. We then evaluate 
these terms by using the eikonal approximation to the Moller operator. Pauli exclusion effects 
are neglected in this preliminary treatment of the FSI. A rough estimate of exclusion effects is 
made from a medium, modified, two-body cross section. Proton projectiles are then compared 
with experiments based on the shell model in a harmonic oscillator basis. Included are the s, p, 
and d shells that allow calculations for targets up to 40 Ca. 

Multiple Inelastic Collision Series 

In the eikonal coupled channels (ECC) model (refs. 11-13), the matrix of scattering ampli- 
tudes for all possible projectile-target transitions is given by 



where barred quantities represent matrices, b is the impact parameter vector, q is the momentum 
transfer vector, and k is the projectile-target relative wave number. In equation (1), Z is an 
ordering operator for the ^-coordinate that is necessary only when noncommuting two-body 



interactions are considered. The eikonal phase elements are defined by matrix elements of 
arbitrary projectile-target states of the following operator: 


*<») - E /” - / A ^ 


( 2 ) 


where a and j label the projectile and target constituents, respectively; r is the internal nuclear 
coordinate; 77 is the projectile-target separation with 77 = (b , 2 ); is the nucleon-nucleon 
(NN) amplitude; and is the NN relative wave number. 

When treating inelastic scattering, we assume that the off-diagonal terms in X (denoted by 
X Q ) are small compared with the diagonal ones, X D \ then we expand f in powers of X Q to 


*■>-£/ 


d 2 b e*q b e**»( b ) 


E {i«£)r} 


(3) 


We also will make the assumption that all the diagonal terms are represented by the ground- 
state elastic phase X. Using equation (3), we sum over target final states X (continuum) to find 
the inclusive angular distribution for the projectile when its mass remains unchanged as in 

= XX fdh dh’ e'«< b “»'> [X(b)-XUb')] 

IN (2^) 2 J 

x E E rV < °P°t\ Rb)] m | 0 pX X X0 P \ [-iX f (b')] m | 0 p 0 T > ( 4 ) 

X/0m=l 

Equation (4) only allows for a study of the momentum transfer spectra of the projectile. However, 
in any consideration of the projectile energy loss, energy conservation must be treated. Based 
on continuum states for the target final state, energy conservation leads to 


da \ 

dn) 


d?a \ k 2 

dQ dEpt J (27t) 2 


f dHdh' e i[^(b>-XU b ')] £ ^(b.b'.u,) 

^ m= 1 


and 


da \ 


dEpf 


) = fdhe- 2l ™ X WJ2w m (b,b,u) 

/IN J m=l 


(5) 

(6) 


where Ept is the energy of the projectile in the final state, a ; is the projectile energy loss, and 
we define 


1 f m 


dkj 

(2tt) 2 


6(E f - Ei) < 0 P 0 T | [X(b)] m |0pkj > 


x < kjOp| X (b ) |0p0p > 


where kj is the wave number vector of a knocked-out target nucleon. 


(7) 
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The inelastic collision series of equation (5) is expected to converge fairly rapidly. In the next 
section we consider the evaluation of this series for an uncorrelated target wave function. 

Collision Terms in Plane Wave Approximation 

We first consider the evaluation of the collision terms W m using plane waves for the final 
continuum states of the target. The projectile motion is treated in the coherent approximation. 
(See ref. 1.) We consider a three-dimensional representation of the collision terms rather than 
just the transverse terms considered in references 2 and 3. The effects of correlations are not 
treated herein. 

The first collision term is written 

A ^ a^t P°° r°° f 

Wi (b, b', a ;) = - P ./ 2 / dz I dz r / dq dq' 

(z7rj /e NN d— oc J — OO J 

x F(q)F(q )f^(q)f^ N (q ) J ~ Ek)Go T k( c l) (8) 

where G 0r k is the transition form factor of the target and A P and A T are the projectile and 
target mass numbers, respectively. It is helpful to change variables as 


a = jfaW) 

(9) 

P = q - q' 

(10) 

x = r — r' 

(11) 

y = \ (r + r') 

(12) 

V 

1 

pr 

II 

(13) 

s = \(v + v') 

(14) 


with the transverse parts denoted and Sj_, respectively. The first collision term is rewritten 
using equations (9)-(14) as 


Wi (Rj_, Sj.,0/) = f dz dz' d 3 a d 3 p e ia R e i/3 S 

(2 TO J 

X A (“ + f ) (" ~ f) 

where we have defined 


>l(q) = F(q)/ NN (q) 

and 

Ri(a,(3,uj) = J - E k )G ok ^ ^ 


(15) 


(16) 

(17) 
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Introducing the Fourier transform pair 


= ( 18 ) 

and 

Rl (a,/3, t) = J duj e~ iu>t Ri (a, /3, w) (19) 

allows us to evaluate the energy- conserving delta function in equation (17). (See ref. 14.) For 
the target nucleons, we assume that 


k 2 

Ek ~ 2^ + ' S > 


( 20 ) 


where mjy is the nucleon mass, is the binding energy, and equation (19) is 


~ . f dk 

= j ^3 dx dy 


-ik z t/2m N e zk-x 


x e* a x e i0 y p{ y + |.y-|) 

where the density matrix is p{ r, r') and is defined by 

p{ r,r') = $(r) 4>V) 


( 21 ) 


( 22 ) 


and $ is the ground-state single-particle wave function. We then find 

Bl ( a , J 9, 0= ^ft/rf» I d s 1/ ^*-‘« i «»j„ ( e 1 x ) p(y+5. y -f)e(a.- tBl ) (23) 

where j 0 is a spherical Bessel function, 0 is the unit step function, and 

£l = y/2rn N (u ~ e Bl ) (24) 

The higher order terms are more difficult to treat because of the enumeration of projectile 
and target intermediate states. A first approximation is to assume that the projectile remains 
in the ground state throughout the collision (coherent projectile approximation). 

Using similar coordinate changes as described above, we find the mth-order collision term as 

*7 fi b ** 

X e iQ J R e^r s Aj (otj + a] (a,- - ?f) 

X Rmi&l ? * • • ? j ■ ■ * i fifTh (^) 
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where 


R m (ai, . . . , a m , 01, ■ ■ ■ , Pm,u) = 


2 


^ \3ro/2 
2n) 


/SI 


cfixj dPyj e* 


*j x j y j 


x p 



f .,l *3m/2— 1 

2J_ | sw 

2 / J / \3m/4-l/2 

N 


x J'lm/2-l 


l2m N (i 


UJ 


€ B, 


.)£ 

j=l 




6 ( w - £ B m ) ( 26 ) 


where i? 7n = 0 for lj < e# m . The solutions for the rath-order terms in equation (26) result from 
the Fourier transform of the temporal response. Because we keep the longitudinal momentum 
transfer in the response (as opposed to the approach of cylindrical geometry in references 2 
and 3), the order of the Bessel function in R m differs from references 2 and 3. For forward- 
peaked wave functions, we approximate 


^m(°q » ■ • • > Pi j • • • > <*0 — — e# m ) m R\ (txji ^ x j') (27) 

j=l ^ v / 


such that 


Wm (Rj_ , S j_ , u;) — 


C m {v ~ tB m ) 
( m !) 2 


m- 1 


wi 


(r_l,s Xj -&L^ 


vWJ 


(28) 


2 

where Ci = 1, C2 = f , C3 — yjj^, and C4 — Equation (27) is found by considering the 
Taylor series for J 3m y 2 -l* The effective energy shift in equation (27) for the m > 2 approximation 
and the coefficients C m differ from references 2 and 3 because of the longitudinal contributions 
to R m that are included here. We then have for the energy loss spectra (eq. (5)) in a coherent 
projectile model, 


and 


dr a 


) = [ Sr SS e i[VRx+Sx/2)-V(R x -S ± /2)] 

An (2?r ) 2 J 


d£l dE pi / j]\r (2 tt) 

A t 


E 

ra=l 


Cm{u ~ e B m ) 


m— 1 


(ra !) 2 




(29) 


da 


dE 


P f 


= J d 2 S e" 2 Im 


At Cm(u-e B J m - 1 


E 

777 — 1 


{ml) 


n2 


w 1 




(30) 
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The coherent approximation assumes that the projectile remains in the ground state through- 
out the scattering. The leading-order correction to the coherent terms occurs in the collision 
term W‘j and corresponds to the following replacement of W 2 from references 2 and 3: 


F ( ai + F (a! - F (a 2 + y ) F (a 2 - f ) 
4 | F( 2ai) + (. A P - 1) F (an + y) ^ («l - y) 
F(2a 2 ) + (A P — l) F (a 2 + y) ** («2 - y)] } 


(31) 


which follows from using closure on the projectile intermediate states. Physically, equation (31) 
allows the projectile to dissociate in the intermediate state. Further modifications are necessary 
when correlation effects are treated. 


Final State Interactions 

The target transition form factors will describe the effects of the FSI between the unobserved 
ejected nucleons and the recoiling target nucleus. We now consider these effects using the 
eikonal form of the optical model. Because the measurement is of the fast primary, we expect 
the details of the FSI on the primary wave function to be small. Therefore, we will introduce 
several approximations to obtain a tractable solution. 

The transition form factor of the target appearing in the first-order response is given by 

Go rkl (q) = < Or I f ' q r I > ( 32 ) 

where is the outgoing scattering state. With the Moller operator the transition form 

factor is written using plane-wave states as (ref. 4) 

G' 0 r k 1 (q)=<OT|e iqr oS c ; ) l k i > ( 33 ) 


The Moller operator is related to the Green function ^ and to the transition operator T as 

6£ )= 1 + 9f ] T (34) 

Using equations (33) and (34), we can separate the first-order response function into three terms 
corresponding to the plane- wave response, elastic distortion in the FSI, and inelastic reaction in 
the FSI (cascade). Thus, 


i?l(q,q'^) = i?r W + A? W + 

The plane- wave term was described above. For the DW term, we have 

< Or | e iq r gi^T |k x > < kj | e - iq ' r ' \0 R > 
+ < Ofl | e iq r | ki > < ki | gi~^T^ e~ iq r \0 R > 


(35) 


r bw 




( 36 ) 
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where |0/j > is the ground-state wave function of the recoil nucleus. The cascade term describes 
a new inelastic collision series of the ejected nucleon with kj reacting on the target recoil given 
by 


R 


f d 3 k 

1 ~hJ (M 


.IN 


* n 


£=2 ' 


3 = 2L 


d?kj 

(2^ 


tlv-E ki Ek j 


3 = 2 


X 


X 


< Or I e^- r 9 { o ] T | kj kj > 

3 = 2 

£ 

< ki n k J I 9 { <T )] T ] e-^' r ' \0 R > 

3 = 2 


(37) 


We next consider the evaluation of formal equations (36) and (37) using the eikonal form of the 
optical model. 


In the optical model (refs. 5-7), the Moller operator is expressed by the matrix 

nt _) = e ***W) 

k i 


(38) 


where s denotes the transverse component of r and the subscript R indicates the coupling phases 
for the recoil system with matrix elements 


X 


(-) _ (At ~ i) 


R 


/ +oc r 

dzjdq c ?q r Go Tkl (q)/ NN (q) 


(39) 


(27r) 2 fc NN 

where the energy dependence of X L '' is determined by kj rather than by the beam energy. 

The diagonal part of ' determines the PW and DW response terms that we combine as 

^ L (q, q', A = -Rf W (q, q', w) + R? W ( q, q', u) (40) 

If we neglect incoherent contributions to the elastic distortion, we have 

d3k i 


0 iq r —iq' r' 


^F L (q,q',w) = J s(u- E kl ) J dr dr ' 

x e _tk i( r_r, )p ( r ,r') expji [x^ -) (s,z) - X^ -)t (s',z') } 


(41) 


Expanding the phase of the distorted wave in equation (41) about r — r' = 0 and keeping 
only the first term were shown to provide accurate approximations for the distorted wave in 
reference 15 and are used here. Thus, 


Rf L ( a , (3,u;) -0^ 6 A-E k )J d 3 x d 3 y e ia x 

x e* k x p (; 


X X 

,y+r y ~ 2 




( 42 ) 


The calculation of X^, ^ is described in the appendix. 


7 



We make one further approximation by noting that the two-body parameters in ^ change 
smoothly with the energy loss and provide the only dependence on kj for this phase. Thus, we 
write 


Rf L (q,q',u;) = ^ J d 3 x d 3 y e ia * joUix)Q {u - e Bl ) 


( XW X \ _2 Im X ( „ _) (y) 

X K y+ 2 )( y " 2 ) e 


(43) 


with the effective energy used for evaluating X^ ^ given by E = u — e^. The Pauli exclusion 

effects should be more important in Xr than the P-T elastic coupling and will be approxi- 
mated by the effective two-body cross section discussed below. 

We next consider the inelastic part of the response function. As in equation (3), we expand 
nj c " ) into diagonal and off-diagonal parts. The off-diagonal terms correspond to inclusive 
reactions between ki and the target recoil. The first term corresponds to ki ejecting a second 
nucleon from the target where we are ignoring low-lying excited states. That term is given by 




<>' 1 "> - / (§T f (“ - - - £ k2 ) / * * * r 


e* k r(r-r') p ( r?r ') e 


X/ ? (S,2)-X^(S , ,2') 


x < Or | Xr(s,z) | k 2 > < k 2 | Xfl(s',z') |0/j > 


(44) 


Using approximations similar to those made in equations (27), (42), and (43), we reduce 
equation (44) to 


- [J (» - <Bj )] ^ /<<* dy e iax e‘»y 

x p ( y + f 1 y ~ f ) i°(& x ) w R ) ( x ’ y- &) 

where we define the collision term between recoil and nucleon knockouts as 

wjfW.fc) = (i4r ^ - f°°dz r dz' [ d 3 a 1 d 3 0 e ia '*W jW) 
'' (2-7r) 4 k NN dx z Jy z J 

x fNN( Q/ + /3V2)fNN( Q/ — /5 / /2)^?i > '^ , 0 , £ 2 ) 


(45) 


(46) 


where x(z) = (xj_,z) and y(z) = (y±,z'). In equation (46) the plane-wave response function 
appears evaluated at the value £ 2 . Higher order terms in the intranuclear cascade could be added 
to the response function in equation (45) in a similar manner. 
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Method of Calculation and Results 

By using shell model wave functions in a harmonic oscillator basis, the density matrix for 
Ap < 40 is found as 


( x 

7+2’ y 


X' 

2> 



„-y 2 !R 2 T e - 


2 /ar 2 


^ a s + 2 Op 




+ a d 



(47) 


where the constant represents the target matter radii and a s , a p , and a,/ are occupation 
probabilities for s p- ■ and (/-shell nucleons, respectively, given by (for .4/- < 40) 


4 

( Ap > 4) 

1 

(A T < 4) 

0 

(Ap < 4) 

Arp 4 

Ap 

S' 

r-H 

VI 

A 

12 

Ap 

(A T > 16) 

0 

(Ap < 16) 

Ap — 16 
Ap 

(Ap > 16) 


(48) 


(49) 


(50) 


In equations (47) (50) we treat the degenerate Id and 2s shells approximately as a single shell, 
denoted d shell, because spin effects are not considered. In table I, values for shell-model 
occupation probabilities are given for several nuclei. 

The plane-wave response function for the s shell is found by using equation (47) in 
equation (23) as 


R s 


rn N a H R T e -R^/3 2 /4 

Cty/ft 


^(« 2 +^ 2 ) S inh(2 R^a^j 


(51) 


Table I. Shell Model Parameters 


Nucleus 

ds 

dp 

dd 

a Rp, fm 

4 He 

1.0 

0 

0 

1.33 

6 Li 

.57 

.143 

0 

2.11 

9 Be 

.444 

.185 

0 

1.79 

12 C 

.333 

.222 

0 

1.69 

16q 

.25 

.25 

0 

1.83 

20 Ne 

.20 

.20 

.033 

2.14 

27 A1 

.148 

.148 

.068 

1.91 

40 Ca 

.1 

.1 

.1 

2.10 


a Values for Rp from references 1 and 16. 
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The plane-wave response function for the p-shell is found by 


m N a p R T f 

Rp -~^r\ 


2- \ R^p 2 + 2 4 (a 2 + £ 2 ) 


sinh ^2Rja£^ 
-4 4 a£ cosh (2 4<*£) | e~ R T e - R r^ 2+ ^ 


( 52 ) 


and the plane-wave response function for the d-shell is found by 

= n^Rj- f I" 7 + 2R 2 (a 2 + ?2) + 2R i (a 2 + ^2,2 

ay/ir I 2 


+ 8 4 « 2 c 2 


+ 


sinh ^2 Rj-ot^j — Rq'P 2 


6 - + R^(a 2 +e) 


sinh 


( 244 ) 


l4’4 sinh ^2 i? 2, — ^Rj + 8-4 a £( a2 + £ 2 )j cosh ^24’ Q £) 

+ 2 4/3 2 «C cosh (2 4 <*) j e-4^ 2 / 4 e -4(“ 2 +^ 2 ) 


(53) 


In an evaluation of the collision terms, the energy loss is taken as Ui — to e b ^ where is 
the separation energy corresponding to the orbit i. (See ref. 10.) Values for eg. for a few nuclei 

are listed in table II as found in references 10 and 16. Analytic forms for and /? IN are not 

possible, so numerical integration is used. 


Table II. Binding Energies of Shell Model Orbits 


Nucleus 

Orbit 

°f Bv MeV 

4 He 

s 

20.5 

6 Li 

s 

26.0 


P 

5.2 

9 Be 

s 

27.2 


P 

18.2 

12 C 

s 

38.7 


V 

17.5 

I6q 

s 

39.0 


p 

18.0 

27 A1 

s 

50.0 


V 

22.0 


d 

15.0 

40 Ca 

s 

51.0 


p 

35.0 


d 

15.0 


“Values for cg 1 from references 10 and 16. 


10 



For Ap < 4 we use F(q) = e R P q2 / 4 i where the constant Rp represents the projectile matter 
radii. The two-body amplitude is parameterized as 


f _ °{P + *)% N -Bq 2 / 2 
/NN " 4n 6 


(54) 


where a is the two-body cross section, B is the slope parameter, and p is the ratio of real to 
imaginary parts of the forward, two-body amplitude. Values for the two-body parameters are 
listed in table III for calculated versus laboratory (lab) energies. 


Table III. NN Amplitude Parameters 


7i ab , MeV 

C Tpp , fill 2 

a np y fm 2 

B ppi fm 2 

B„ p , fm 2 

Ppp 

Pnp 

558 

3.85 

3.58 

0.12 

0.12 

0.40 

-0.16 

800 

4.67 

3.78 

.16 

.16 

.18 

— .33 


In figures 1-3 we show the diagonal part (/3 = 0) of the response function for 4 He, 12 C, and 
27 A1 at several values of a versus the energy loss. Binding energies assumed for these nuclei are 
listed in table II. The multiple shell structures for 12 C and 27 A1 are seen at lower values of a in 
figures 2 and 3. In figures 4 6 the off-diagonal part (/? ^ 0) is displayed for a = 1.5 fm -1 . The 
breaking of translational invariance in finite nuclei leads to a nonzero response for the off-diagonal 
components. Figures 4-6 show that the off-diagonal terms quickly dampen as we increase the 
target mass number, and the results for nuclear matter should be accurate for nuclei heavier 
than those considered here. Calculations that include ground-state correlations by Alberico et al. 
(ref. 17) also suggest that damping saturates at about Aj = 40. Figure 7 illustrates the complex 
structure of 27 A1 caused by the multiple shell structure and the off-diagonal response. 

The response terms with contributions, for the FSI (eqs. (43) and (45)) contain the free two- 
body amplitude evaluated at an effective energy E that was determined by the projectile energy 
loss as E = w — . For the relatively low values of E, we should expect significant deviations 

from the impulse approximation. We follow the usual approach of replacing the free two-body 
amplitude in equation (54) with an effective one that approximately treats some of the medium 
modifications. Smith and Bozoian (ref. 18) provide the parameterization 

<Ti£i64Mev) 

_ ^ 0.6 T l (T l > 164 MeV) 

where Pi and Ti are the nucleon kinetic momentum and energy, respectively, in the laboratory 
that correspond here to E. The slope parameter is assumed as zero, wdiich corresponds to the 
isotropic scattering that is approximately true at low energies. Equation (55) is based on optical 
potential studies. In reality, & should contain a dependence on the target density and will differ 
for elastic or inelastic collision terms because the number of off-the-mass shell particles is not 
the same in these two cases. 

Comparisons of the distorted- wave and inelastic response corrections with plane-wave re- 
sponse calculations are shown in figures 8 11 for 12 C. The DW response (dashed line) is lower 
than the PW response (solid line) at small values o.f a. The shapes of the response functions for 
the diagonal terms show only a slight shift in the position of the quasi-elastic peak because of 
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the inelastic term. In contrast, the off-diagonal contributions (figs. 10 and 11) are substantially 
modified. Here, the distorted waves alter the interference patterns because of the multiple shell 
structure. Calculations for heavier targets were not performed at this time because evaluations 
of inelastic response function are extremely complex for the higher shells. 

Calculations of momentum spectra are shown in figures 12-15 at several angles along with 
comparisons with experimental data (ref. 19) for inclusive proton scattering on ^Li, ^ 2 C, 2 7 A 1 , 
and 40 Ca at a beam energy of 800 MeV. At the lowest momentum values, the effects of pion 
production are seen but are not included in our calculations. The calculations in figures 12 15 
are made with the plane- wave response model and include contributions up to the fourth order. 
The second- and higher-order terms are small; however, these terms become more important with 
increasing target mass and momentum transfer. In figures 16 and 17, a similar comparison is 
made with the experimental data of reference 20 for protons at 558 MeV. Clearly, the dominant 
contribution for nucleon-nucleus collisions is the first collision term, which suggests that an 
improved description of the data will require going beyond the independent particle model used 
here to include correlation effects. This dominance of the first collision term contrasts with our 
previous results for composite projectiles (refs. 2 and 3) in which coherence effects are important 
for multiple scattering terms. 

In figures 18 and 19, calculations are shown for 12 C and 27 A1 that include the effects of final 
state interactions in the first-order response functions for 800-MeV proton beams. The inelastic 
collision term is seen to contribute to the distribution in a manner close to the plane-wave 
double collision term; this contribution is expected because the terms have the same dependence 
on energy loss and ground-state wave function. The combined effects of the distorted wave, 
ejectile, and recoil reduce the peak of the cross section with a slight shift in position. The 
inelastic collision term contributes primarily in the dip region between the quasi-elastic and 
pion production peaks. 

Concluding Remarks 

A formalism for describing the energy loss spectrum of fast ions in nuclear collisions at high 
energy was developed that includes the effects of final state interactions. Calculations in the 
independent particle model with harmonic oscillator wave functions suggest that a signature of 
the intranuclear cascade is seen in the dip region between the quasi-elastic and pion production 
peaks. The effects of the nuclear medium and the Pauli blocking are expected to be important in 
providing a more complete treatment of cascade effects than discussed in this paper. However, 
the present results indicate that cascade effects can be treated in an approach that is similar to 
the quantum mechanics used here. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
June 14, 1993 
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Appendix 
Calculation of X^ R ^ 

The outgoing distorted phase for the elastic coupling of the ejected nucleon with the core is 
written 

x (-) (v) _ (At ~ 1) 

* ( } ( 27 r) 2 kN N J vt 

The form factor is found from the one-body density matrix in equation (47) as 


2 ' f dz f d3q eiqT C(q)fNN(q) 
2 k;™ Jv? J 


where 


G(q) = (Ci + C 2 q 2 + C 3 q 4 ) e~ R W / 4 


C\ — a s + 3 a p + 6 a-d = 1 


c 2 = -m 


^ - a p + 3 dft'j 


Ri 

c 3 = -f a r 


Using the two-body amplitude of equation (54), we have for the imaginary part of Xr 


(-) 


where 


Ri B 

^X +2 


Integration of equation (A5) leads to 


Im X«->(y) = (At ~l )<r e'd/"* 


Air ( \ y 2 zt 4^ 


(! 4 )- 2 -(l 4 ) 


where T(a, v) is the incomplete gamma function defined by 

roo 

r(a,ti) = / t a ~ l dt 
Ju 

and the coefficients A\,A 2 , and A 3 are defined as 

4 _r , C2/3 y±\ C 3 f 15 5 yj yj ' 

1 1 w (^2 4 w J re 2 y 4 4 in 16tn 2 


(Al) 


(A 2 ) 

(A3) 

(A4) 

(A5) 


im xr){y) = r dz i d3q eiqy ( Ci + ° 2q2 + c3 ' 4 ) e-w92 (A6) 


(A7) 


(A 8 ) 


(A9) 


(A 10 ) 


4 w 2 \ 4 8 w 


(All) 
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and 


(A12) 


4 C 3 

3 16w 4 

The energy variation of Im \ jy ' is through the two-body parameters cr and B that are isospin 
averaged and evaluated at the energy w - in equation (43). 

The coupling between the projectile and target can be evaluated from equation (A8) by 
replacing ( A T - 1) with A P A T , letting y z -> -oc, using two-body parameters appropriate for 
the beam energy, and redefining w as 

I$p Up B / * i q\ 

w = -r + -r + it ( A13 ) 

4 4 2 

where we assume a light projectile when a Gaussian form factor is sufficient. 
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Figure 1. Diagonal part of first-order response function versus energy loss for 4 He with (3 = 0 . 
All dimensions are in fm . 



Figure 2. Diagonal part of first-order response function versus energy loss for 12 C with (3 — 0 . 
All dimensions are in fm -1 . 
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Figure 3. Diagonal part of first-order response function versus energy loss for 27 A1 with /? — 0. 
All dimensions are in fm -1 . 



Figure 4. Off-diagonal part of first-order response function versus energy loss for 4 He at 
a = 1.5 fm _1 and several values of (3. All dimensions are in fm 
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Figure 5. Off-diagonal part of first-order response function versus energy loss for 12 C at 
a = 1.5 fm -1 and several values of /3. All dimensions are in fm _1 . 



Figure 6. Off-diagonal part of first-order response function versus energy loss for 27 A1 at 
a = 1.5 fm -1 and several values of p. All dimensions are in fm -1 . 
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Figure 7. Off-diagonal part of first-order response function versus energy loss for 27 A1 with 
contributions from individual shells, a = (3 = 1 fm“*. 



Figure 8. Diagonal part of first-order response function for 
a — 1.5 fm -1 ; (3 = 0 fm -1 . 


12 C with inelastic corrections. 
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<jl), MeV 

Figure 9. Diagonal part of first-order response function for larger momentum transfer with 
inelastic corrections, a — 2.5 fm -1 ; (3 — 0 fm 1 . 



co, MeV 

Figure 10. Off-diagonal part of first-order response function with inelastic corrections, 
a = 0.75 fm" 1 ; (3 = 1 fm" 1 . 
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Figure 
a ~ 



11. Off-diagonal part of first-order response function with inelastic corrections. 
1.5 fm -1 ; (3 = 1 fm™ 1 . 
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Figure 12. Calculations of inclusive p scattering for ®Li at 800 MeV compared with data of 
reference 19. p + 6 Li -> p + X at 800 MeV. 
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Figure 14. Calculations of inclusive p scatterin 
reference 19. p + 27 A1 — + p + X at 800 MeV. 
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Figure 16. Calculations of inclusive p scattering for 12 C at 558 MeV compared with data of 
reference 20. p + 12 C — > p + X at 558 MeV. 
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Figure 17. Calculations of inclusive p scattering for 27 A1 at 558 MeV compared with data of 
reference 20. p 4- 27 A1 — ► p + X at 558 MeV. 




(b) 0 = 15°. 

for 12 C at 800 MeV that include final state 
p + 12c p + X at 800 MeV. 


(a) © = 11°. 

Figure 18. Calculations of inclusive p scattering 
interactions. Experimental data from reference 19. 




(a) 0 = 11°. (b) © = 15°. 

Figure 19. Calculations of inclusive p scattering for 27 A1 at 800 MeV that include final state 
interactions. Experimental data from reference 19. p + 12 C -*• p + X at 800 MeV. 
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